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Motivated by experiments on non-magnetic triangular-lattice Mott insulators, we study one can- 
didate paramagnetic phase, the columnar dimer (or valence-bond) phase. We apply variants of 
the bond-operator theory to a dimerized and spatially anisotropic spin-1/2 Heisenberg model and 
determine its zero-temperature phase diagram and the spectrum of elementary triplet excitations 
(triplons). Depending on model parameters, we find that the minimum of the triplon energy is 
located at either a commensurate or an incommensurate wavevector. Condensation of triplons at 
this commensurate-incommensurate transition defines a quantum Lifshitz point, with effective di- 
mensional reduction which possibly leads to non-trivial paramagnetic (e.g. spin-liquid) states near 
the closing of the triplet gap. We also discuss the two-particle decay of high-energy triplons, and 
we comment on the relevance of our results for the organic Mott insulator EtMe3P[Pd(dmit)2]2. 

PACS numbers: 75.10.Jm, 75.10.Kt, 75.30.Kz, 75.50.Ee 



I. INTRODUCTION 

In the search for novel phases, frustrated quantum an- 
tiferromagnets have attracted enormous attention. Here, 
the combined effect of geometric frustration and quantum 
fluctuations tends to destabilize conventional magnetic 
order in favor of quantum paramagnetic ground states, 
such as valence-bond solids (VBS) with broken lattice 
symmetries or featureless spin liquids-'^^ 

A prominent example of a frustrated quantum mag- 
net is the spin-1/2 Heisenberg model on the triangular 
lattice. Although it wa.s initially proposed by AndersorP 
that its ground state could be a resonating-valence-bond 
(RVB) spin liquid, it is now established that the model 
displays semiclassical non-collinear 120° order.l^ How- 
ever, modifications beyond the simple nearest-neighbor 
exchange are believed to induce non-magnetic phases on 
the triangular lattice. For instance, it has been arguecP 
that a combination of longer-range and multiple-spin- 
exchange interactions is important near the Mott transi- 
tion to a metallic state, where they induce a spin-liquid 
state as observed in numerical studies of the triangular- 
lattice Hubbard model. Models with spatially anisotropic 
exchange interactions have been studied as well,'^^^^ in- 
terpolating from the triangular lattice to both the square- 
lattice and the decoupled-chain limits, and both VBS and 
spin-liquid phases have been proposed to occur. 

Experimentally, triangular-lattice Heisenberg mod- 
els are realized in a variety of materials, such as 
the inorganic Mott insulators CS2CUCI4 and Cs2CuBr4 
(Ref. il8j and fT9|) and the organic compounds k- 
(ET)2Cu2(CN)3 (Ref. [201 and X[Pd(dmit)2]2 (Refs. [H]- 
[551 with X^EtMesSb or EtMeaP), all showing some de- 
gree of spatial anisotropy. In the latter two, the crystal 
lattice is rather soft: pressure can be used to induce an 
insulator-to-metal transition, with superconductivity ap- 



pearing at low temperature. Spin-lattice interactions are 
also relevant in the insulator, as they potentially relieve 
frustration via dimerization. In fact, such a dimerized 
VBS phase has been reported for EtMe3P[Pd(dmit)2]2, 
likely stabilized by magneto-elastic coupling. 

In this paper, we study the columnar dimer phase of 
a generic triangular-lattice antiferromagnet (AF).'2fil In 
the absence of detailed knowledge about microscopies, 
we focus on a simple, yet very rich, model (Fig. [T^) 
with anisotropic nearest-neighbor couplings and explicit 
dimerization, the latter possibly arising from magneto- 
elastic effects. A determination of the full phase diagram 
of this model (see Fig.[2]for a sketch) is a hard task. Here 
we stick to the more modest goal of characterizing the 
dimer phase and its triplon excitations and locating the 
line of quantum phase transitions (QPTs) to a magneti- 
cally ordered state, where the triplon gap closes. To this 
end, we employ the bond-operator technique introduced 
by Sachdev and Bhatt.'^S It turns out to be important 
to go beyond the linearized (i.e. non-interacting) bo- 
son problem: We discuss various approximation schemes, 
also providing a guideline as to which non-linear effects 
are important in the presence of frustration. 

Among our results is the existence of a commensurate- 
incommensurate transition (CIT) in the excitation spec- 
trum, i.e., the minimum of the triplon dispersion is locked 
at wavevector (0, 0) (in the coordinates of the dimerized 
lattice. Fig. [T]) in some region of parameter space, while 
it moves to incommensurate momenta elsewhere. If the 
triplon gap closes at the CIT line, the excitations be- 
come anom alously soft, leading to a quantum Lifshitz 
point .I^^^^H We briefly discuss the unusual quantum crit- 
ical regime of such a putative magnetic QPT, but also 
speculate about an intervening non-trivial paramagnetic 
phase. Finally, we discuss two-particle (as oppo sed to 
three-particle) decay of magnetic excitations,'2SEIl which 



FIG. 1. (Color online) (a) Lattice structure of the model ([T]), with exchange couplings J, J', and J" shown as black (thin), 
red (thick), and dashed black lines, respectively. Solid (open) circles represent spin (S^) of each dimer. ai and a2 are 
the primitive vectors of the dimerized lattice, (b) Dimerized ground state in the limit of large J'. Here, ellipses represents a 
singlet, (c) Brillouin zone of the dimerized lattice (solid), together with the (hexagonal) Brillouin zone of the original triangular 
lattice (dashed). The open circles indicate the minimum position of Wq [Eq. (17l] for parameters J' = 1.5 and J" = 1. Here 
K = (47r/3a, 0), K' = (27r/3a, 2-K/\/ia), X = (tt/h, 0), Y = (0, 7r/\/3a), and M = (vr/a, 7r/\/3a) with a being the lattice spacing 
of the triangular lattice. 



is generically possible in low-symmetry dimerized mag- 
nets. 



A. Model 

We consider a triangular-lattice Heisenberg model, 

7^ = ^J,,S,.S„ (1) 

<i,i> 

where Si represents a spin-1/2 at site i, and the nearest- 
neighbor exchange couplings Jij equal J, J' , and J" ac- 
cording to the pattern shown in Fig. Hereafter we 
employ J — \. 



B. Phase diagram 

To appreciate the richness of the model Q, it is use- 
ful to discuss a few limiting cases, with their properties 
summarized in the schematic phase diagram in Fig. |2] 

(i) Line J" = 1. On this line, for J' = 1, we recover 
the isotropic triangular-lattice AF Heisenberg model. Its 
ground state has non-coUinear 120° Neel order, with 
commensurate ordering wavevectors K = (47r/3a, 0) and 
K' = (27r/3a,27r/\/3a).^Non-collinear order can be ex- 
pected to extend to small values of J' . On the other 
hand, the limit J' S> 1 corresponds to a paramagnet of 
weakly coupled dimers,'^such that spins which are pair- 
wise coupled by J' combine into a columnar arrangement 
of singlets, as illustrated in Fig. [TJj. 

(ii) Line J" = 0. Model ^ is now topologically equiv- 
alent to the staggered dimerized Heisenberg model on a 
square lattice, recently discussed in Refs. I33landl34l Ac- 
cording to quantum Monte Carlo (QMC) calculations, a 
order-disorder QPT takes place at J'^ = 2.5196,1221 with 
coUinear Neel order setting in for J' < ,][.. The line 



J" = also includes the square and honeycomb lattice 
AF at J' = 1 and 0, respectively. 

(iii) Line J' = 1. Here, Eq. ([T]) is nothing else but 
the spatially anisotropic triangular-lattice AF addressed 
in Refs.i7til7Jn particular, coupled-cluster calculation^iSl 
indicate that the system displays coUinear Neel order for 
J" < 0.8 and non-coUinear spiral order with an incom- 
mensurate ordering wavevector (except at J" — 1) for 
0.8 < J" < 1.8. For J" > 1 the system consists of 
weakly coupled chains, and a coUinear AF stat^^ has 
been a rgued to arise as a result of order-from-disorder 
physics .D^Ji^ However, we note that both non-collineaiti^l 
and disordered (i.e. spin-liquid) ground states^'' have 
also been proposed in this regime. For 0.7 < J" < 0.9 the 
physics is under debate as well: series-expansion studies 
favored a spontaneously dimerized VBS in this region,!^ 
which was not found in the coupled-cluster study.l^l 

In total, the phase diagram of the model ([T]), Fig. [2] dis- 
plays a gapped paramagnetic phase, both coUinear and 
non-coUinear long-range order (LRO), as well as puta- 
tive non-trivial spin-liquid regimes which may or may not 
be adiabatically connected to the one-dimensional limit 
J" oo. 

In this paper, our focus is on the properties of the 
dimerized paramagnet which is adiabatically connected 
to the limit of large J'. As stated above, in a real sys- 
tem, such as EtMe3P[Pd(dmit)2]2 the dimerization may 
arise spontaneously due to longer-range or ring-exchange 
couplings, and will then be stabilized by magneto-elastic 
effects. Leaving a detailed study of the latter for the fu- 
ture, we choose to work with explicit dimerization as in 
Fig. [ipSI Most of our calculations are restricted to the 
parameter range J' > J" and < J" < 1.5. A quantita- 
tive phase diagram obtained from various bond-operator 
approximations is in Fig. [3] and will be discussed in detail 
below. 
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FIG. 2. (Color online) Schematic phase diagram of the 
Heisenberg model ([T]), as function of the couplings J' and 
J", keeping J — 1. For large J', a gapped dimer phase is 
realized (shaded). Various other limits are indicated in the 
figure, see text for details. The multicritical point, where the 
coUinear LRO, non-coUinear LRO, and gapped dimer phases 
meet, is the quantum Lifshitz point. [The classical limit of 
([T]) is discussed in Appendix [A|] 



C. Outline 

The remainder of our paper is organized as follows: 
In Sec. [nj we briefly summarize the bond-operator for- 
malism that is employed to describe the paramagnetic 
phase of the model ([!]) and derive an effective Hamil- 
tonian of interacting triplet excitations (triplons). In 
Sec. Ill we analyze the non-interacting (i.e. harmonic) 
part of this Hamiltonian. We determine the magnitude 
and momentum-space location of the minimum energy 
gap of the triplons. Corrections arising from interactions 
- both cubic and quartic - among the triplons are ana- 
lyzed in some detail in Sec. |IV[ For all levels of approxi- 
mation, the closing of the triplon gap allows to construct 
the phase boundary between the singlet and long-range- 
ordered phases, shown in Fig.|3] In Sec.jVj we discuss var- 
ious aspects of our results, such as the two-particle decay 
of triplons due to cubic interactions, the commensurate- 
incommensurate transition, and the associated quantum 
Lifshitz point. We also comment on some implications for 
the organic Mott insulator EtMe3P[Pd(dmit)2]2- Con- 
cluding remarks close the paper. A summary of the clas- 
sical phase diagram of ([T]) as well as details of the calcu- 
lations are relegated to the appendices. 



II. BOND OPERATORS AND TRIPLON 
EXCITATIONS 

For J' > J, J" it is useful to view the spin pairs cou- 
pled by J' as building blocks of the model. The four 
states per such dimer can be efficiently represented using 
bond operators.'^ which naturally lead to a description 
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FIG. 3. (Color online) Quantitative phase diagram for the 
Heisenberg model |T]), obtained using the bond-operator for- 
malism. The dimer phase boundary (□) has been obtained 
from the HF-cubic approximation described in the paper; for 
comparison we also show the results of the harmonic mean- 
field^^©) and harmonic spin-wave approximations (dashed 
blue) .'^ The dashed-dotted line separates the commensurate 
(C) and incommensurate (INC) regions of the dimer phase 
at large J'. It continues into the dotted line which separates 
the phases with collinear and non-collinear LRO at small J' 
- note that this line is an estimate based on series-expansiorP 
and classical-limit results (Appendix [A]) Also shown are the 
isotropic point J' = J" = 1 (a) and the point J' — 1, 
J" = 1.05 ( ) for EtMe3P[Pd(dmit)2]2.'^(Here, longer-range 
or multiple-spin interactions can be expected to be relevant in 
addition.) Finally, on the J" — axis we have indicated the 
QMC value of (♦) for the square-lattice staggered dimer- 
ized Heisenberg model.'^ 



of the elementary excitations of the paramagnetic phase 
in terms of bosonic spin-1 modes, dubbed "triplons". 



A. Bond-operator representation 

To introduce bond operators, the triangular lattice of 
spins is re-interpreted as rectangular lattice of dimers. 
Fig. [T]d, with sites i. For each dimer, one introduces 

bosonic operators {s1,4q} ("^ — •^)' which create the 
dimer states out of a fictitious vacuum. Explicitly (and 
omitting the site index i), \s) = s^O), \a) = <jj|0), where 
k> = (in> - \it)_)/V2, \x) = {- |n> + \ii))/V2, \y} = 
^(Itt) + \U))/V2, \z) = (in) +\it))/V2. The Hilbcrt- 
space dimension is conserved by imposing the constraint 



(2) 
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on every site i. The original spin operators and of 
each dimer are given by 



(3) 



where eap-y is the antisymmetric tensor with e^yz — 1 and 
summation convention over repeated indices is imphcd. 



B. Effective theory for triplet excitations 

Expressing the Hamiltonian ([T]) in terms of the dimer 
spins (S^ and Sf) yields 

n = J2 [J'^l ■ + J" (s^ • si, + s2 . s^i)] 

i 

+ JJ2SISI^. (4) 

Here n = 1,2,3 corresponds to the nearest-neighbor vec- 
tors 



Tl 



ax, 



T2 



ay, 



T3 



^a{x + V3y), (5) 



with a being the lattice spacing of the original triangular 
lattice (in the following a = 1). Substituting Eq. ^ into 
Q yields a Hamiltonian of the form 



(6) 



where Hn contains n triplet operators (see Appendix |b] 
for details). 

In the paramagnetic phase realized for large J', the sin- 
glet can be viewed as "condensed" ^ formally s ■ , Si — ?> 
VIVo in Eq. As a consequence, one ends up with 

an effective Hamiltonian solely in terms of the (bosonic) 
triplet operators t^^. In the spirit of mean- field the- 
ory, the constraint ^ is implemented on average via 
a Lagrange multiplier ^. Both fj, and iVo will be self- 
consistently determined. One expects that A'q ~ 1 in the 
limit J' ^ J, J" . (Other bond-operator schemes will be 
discussed in Sec. II C| below.) 

After performing a Fourier transform, i.e., = 
A^'^-^/^Vj^ exp(— ik • Rt)*ka' terms of the Hamilto- 
nian ([6]) read 



En = 



-3J'N/8- fiN{No- l)/2, 

Ma^ka+^Bk(iL*L+H.C.) 

k L 

^-^Ea/JA Ck-p tl_pjlptwx + H.C., 
p,k 

1 V -v I I 



2N' 



(7) 
(8) 
(9) 



q,p,k 



Here A''' = N/2, with A^ being the number of sites of 
the original triangular lattice, and the momentum sum 



runs over the dimerized (rectangular) Brillouin zone. The 
coefficients Ak, i?k, Ck, and 7k are given by 

J' 



Ak = — - ^i + Bu 



(2 J" — 1) cos kx — cos{V^ky) 
— cos{kx + VSky) 



/No 



iinkx + sm{V3ky) + sin(fcj; + Vsky) 



7k = 



(2 J" + 1) cos k^ + cos{V3ky) 
+ cos(A;2, + y/Sky 



(10) 

(11) 

)] , 
(12) 

(13) 



A few remarks about the general structure of the ef- 
fective Hamiltonian ^ are here in order: The bond- 
operator approach has some parallels to the Holstein- 
Primakoff approach to ordered magnets, with the differ- 
ence that one considers fluctuations above a quantum 
paramagnetic state instead of a Neel state. A cubic 
triplet term ([s]) is present in many low-symmetry situ- 
ations, including the model considered here and also the 
square-lattice staggered dimerized Heisenberg model of 
Ref s . [551 and IMl This is to be contrasted with spin waves 
where a cubic interaction term - describing two-particle 
decay of transverse magnons - is only present for non- 
coUinear spin structures while it vanishes for collinear 
ones (see Ref. |3S]for a discussion). 



C. Alternative bond-operator schemes 

The procedure discussed in the previous section, where 
the constraint ([2| is included into the description via 
a Lagrange multiplier and Nq is self-consistently deter- 
mined, is the one originally proposed by Sachdev and 
Bhatt.!^ However, alternative methods to deal with the 
constraint ([2| can be found in the literature. Let us 
briefly summarize two of them. 

Chubukov and Mori'2Sl resolve the constraint via s — 
(1 - Xtitay/^ where A is an artificial control parame- 
ter, with A = 1 corresponding to the physical case. The 
square root can now be expanded to obtain a Taylor se- 
ries in the control parameter A - a procedure similar to 
the 1 /S expansion in conventional spin- wave theory. This 
generates a Hamiltonian with triplon terms up to arbi- 
trary order, which can be analyzed order by order in A. 
Finally, one arrives at physical results by taking the limit 
A^ 1. 

Kotov et alW^ instead implement the hard-core con- 
straint explicitly in the following way. First, the t opera- 
tors are re-interpreted as creation operators of triplons on 
top of a singlet background. This is formally equivalent 



to setting s. 



1. The resulting Hamiltonian con- 



tains terms up to quartic order; its harmonic piece, 7^2, 
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FIG. 4. (Color online) Upper row: Triplon dispersion relation at the harmonic level (solid black) along paths in the dimerized 
Brillouin zone (Fig.mc), for (a) J' = 3.0, J" = and (b) J' = 1.5, J" = 1. For comparison, the dispersion obtained within 
the HF approximation for the quartic term (thin green line, see Sec. IV A I is also shown, indicating that the corrections from 
the quartic term are minor. Finally, the dashed red line represents the bottom of the two-particle continuum at the harmonic 
level. Lower row: Contour plot of the triplon energy at the harmonic level, for (c) J' = 3.0, J" = and (d) J' = 1.5, J" = 1. 



is analogous to linear spin-wave theory. The constraint 
is converted into the inequality tla^ia 5: Ij which can 
be imposed using an infinite on-site triplet repulsion term 



2N' 



■u 



E 



U ^ oo. 



Hjj, which provides the main contribution to the renor- 
malization of the noninteracting triplet energy, is treated 
using a self-consistent ladder summation of diagrams 
(Brueckner approximation) . 

Both schemes were applied to the non-frustrated 
square-lattice bilayer Heisenberg model, with the Brueck- 
ner approach'^ giving very accurate results, e.g., for 
the location of the phase boundary. Results of similar 
quality where obtained for other models with coUinear 
spin correlations.'^ In contrast, we have found that the 
Brueckner approach is less well suited for the present 
triangular-lattice model: Following Ref. [371 we have im- 
plemented the Brueckner approximation for the hard- 
core triplon repulsion and included a non-self-consistent 
treatment of the cubic term to second order. In the 
resulting phase diagram, the stability of the paramag- 



netic phase is clearly overestimated, i.e., we find the 
gapped paramagnet to be stable even at the isotropic 
point J' = J" = 1, with the critical w 0.75 at J" = 1. 
We suspect that the combined effect of the hard-core and 
cubic terms in the presence of non-collinear correlations 
requires a more accurate approximation, but have not 
explored this route further. Therefore we resort to the 
mean-field-based approach of Ref. [17] whose results are 
presented in the body of the paper. 



III. HARMONIC APPROXIMATION 



The lowest-order approximation to the triplon dynam- 
ics consists in keeping the quadratic term 7^2 of the 
Hamiltonian only, which describes the physics of nonin- 
teracting bosons. V.2 can be diagonalized with the help 
of the Bogoliubov transformation 



(14) 



6 



(a) 



o -0.2 



-0.4 



-0.6 



1.7 



harmonic 

O HF 

O HF-cubic 



0.6 



O 



O /• o 



n l o-9-<>«^ 





1 1 1 1 1 1 1 1 

O-O-O0O 


1 1 1 1 1 1 










°o ^ 






, 1 . , , , 1 , 


- * * ^ ■ 

. 1 , ®-~-e-i 



0.7 



0.8 
J" 



0.9 



(b) 



J' = 3.0 



harmonic 

O HF 
2 1- O HF-cubic 



o 



1 : 
Otj 



9 




Oil- ■ 

-0.2 - 

o 

a'^-0.4r 
-0.6 r 
-0.8^ 








0.5 



6- 4) 



1.5 



FIG. 5. (Color online) x and y components of the momentum Qo (see text for definition) as a function of J" at the harmonic 
(dashed), HF (O), and HF-cubic (0) approximations for the lines J' — 1.7 (a) and 3.0 (b). Non-zero values of Qo correspond 
to incommensurate correlations. 



One finds 



where 



(15) 



k a 



Eo = ^p'NN„ ~ l^iN{No - 1) -{- 1 ^ (wk - Ak) (16) 



is the ground state energy, 



(17) 



is the energy of the triplet excitations, and the Bogoli- 
ubov coefficients in Eq. ([T4| obey 



Wk, -^k = ±1/2 + ^k/2t^k and u^Vk = Sk/2wk- (18) 
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FIG. 6. (Color online) x and y components of the momentum 
Qo (see text for definition) as a function of J' at the harmonic 
(dashed), HF (O), and HF-cubic (0) approximations for the 
line J" = 0.8. In the strong-dimer limit of large J', the 
harmonic Qo is approached for both the HF and HF-cubic 
results. 



From the saddle points conditions dEg/dNo = and 
dEo/dfj, = 0, self-consistent equations for /i, and Nq follow 
(see Appendix [C] for details) 



3J^ 



3 

No 



Nn 



^ BkWk (wk - Wk) , (19) 

k 

-i-]|7E-^ = i-^E(^Lw. (20) 



k a 



Once II and Nq are known, the triplon energy (17) is 
completely determined. 

In Fig. |4j we show the triplon dispersion relation for 
two sets of model parameters inside the disordered phase. 
For the case J' — 3 and J" — 0, the minimum gap is lo- 
cated at r, the center of the dimerized Brillouin zone, 
while for J' = 1.5 and J" = 1.0, the minimum gap is 
at incommensurate momenta, Qo = ±(27r/3, — tt/S-s/S)- 
Indeed, within the harmonic approximation, the momen- 
tum Qo associated with the triplon gap can be analyti- 
cally calculated. From the solution of VqCJq = VqSq = 
0, one finds that it depends only on the coupling J", 
namely 



Qo 



with 



(0,0), 

±(^xm7 Qyni) 



2 arccos 



< J" < 0.75, 
J" > 0.75, 



(21) 



4J" 



= arccos _ 2 



(22) 



Note that Qo continuously moves from a commensu- 
rate, (0, 0), to an incommensurate point, ±(gxin, 9ym), as 
J" increases, defining a commensurate-incommensurate 
transition (CIT) in the excitation spectrum at J^ij^ = 
0.75, see Fig. [5j At the present harmonic level, Qo only 
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FIG. 7. (Color online) Parameter A^o as a function of J' for 
the lines J" — and J" = 1 at the harmonic and Hartree- 
Fock levels. 



depends on J", but this does not hold once corrections 
are included, see Se c.|IV[ Further aspects of the CIT will 
be discussed in Sec. IV bI below. 

Parenthetically, we note that - at the harmonic level - 
one recovers the spin- wave-type approximation of Kotov 
et alW^ (see Sec. II C ) if, instead of calculating fi and Nq 
self-consistently, one sets fi — — 3J'/4 and Nq — 1. 



IV. TRIPLON INTERACTIONS 



Both the cubic and quartic interaction terms, ([8| 
and 7^4 ([9]), renormalize the triplon energy and, con- 
sequently, shift the phase boundary of the paramag- 
netic phase. Moreover, the cubic term also induces two- 
particle decay of triplons as discussed in Refs. and 

EH 

In order to include triplon-triplon interactions into 
our description, we treat the effect of 7^4 at the mean- 



field level (Sec. IV A I and the one of "Hs perturbatively 
(Sec.[rVBl). 



Here 



kq 



TV' ^ 
q 



^Arcos(k.r„), 

n=l 



(24) 



N' ^ 
q 



7q-kl 



3 

E 

n=l 



i?«''cos(k.x„), 



with 7k being the bare quartic vertex ( |13[ ), Uq and Uq 
the Bogoliubov coefficients (see definition below) , t„ the 
nearest-neighbor vectors ([5]), and the coefficients A, 
read 



HF 



Af' = (J" + 2) 



k 



cos 



^ TV' ^ 

k 



(\/3/cj,)? 



k' 



(25) 



Similar expressions hold for but with it^ 
The final Hamiltonian is now given by 'H. 



Eq 



Mk^k- 



H2 + "HJ^. It is quadratic in triplet operators and can 
also be diagonalized by the Bogoliubov transformation 
( 14 1, namely 



n = E, 



HP 



E 



(26) 



In the above expression, the renormalized triplon energy 
Wk is equal to Eq. ( 17 1 but now ^k — > ^k = ^k+^k^ and 
Bk Bk 



Bu 



. The ground state energy Eq 



HF 



is similar to Eq. ( 16 ) apart from the replacements ojk ^• 
Wk and Ak — >■ Ak and the inclusion of E]^^ . Finally, 
the Bogoliubov coefficients and the analog of the self- 
consistent equations (19 1 and (20) now read 



-2 -2 
Mk,t^k 



±l/2 + ik/2wk, SkWk = Sk/2wk, (27) 



^J' = — ^ + ]^ E ^k«k (wk - Uk) , (28) 



A. Hartree-Fock approximation 



No = i-—yv 



(29) 



We treat the quartic triplon interaction within the 
self-consistent Hartree-Fock (HF) approximation. It is 
straightforward to show that Eq. (|9l assumes the form 
(for details see Appendix A of Ref.l35|) 



HF 



'-ka'-kQ 



H.c. 



(23) 



Note that _Bk (not Bk) enters the momentum summation 
in Eq. (281 (for details, see Appendix [C]) . In addition to 



ji and A'o, now 3 and -81^2 3 are _also self-consistently 



calculated. The set of equations (23)-(29) constitutes 
HF approximation. 

The resulting triplon dispersion is included in Fig. [4j 
which shows that corrections arising from "H^^ are small, 
except for a slight upward renormalization of the gap. 
For a fixed J', the minimum wavevector Qq displays the 
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(a) 



(1) 



(2) 






(e) ^ 



FIG. 8. (a) Renormalized cubic vertices ri(k, p), (1), and 
r2(k,p), (2). Here (a,/3,7) = {x,y,z), {z,x,y), and {y,z,x). 
(b)-(e) Lowest-order diagrams resulting from the combination 
of vertices (1) and (2) which contribute to the normal self- 
energy E3(k,a;). 



same qualitative behavior in terms of J" as found in the 
harmonic approximation, see Fig.[5j In the incommensu- 
rate regime, Qo numericahy deviates from the harmonic 
result for small J' (but recovers the harmonic results in 
the large- J' limit), see Fig.[6j Moreover, the line marking 
the CIT is slightly shifted as well, see Fig. [Sj 

Fig. [t] shows the evolution of A'o as a function of J' for 
fixed J", again comparing the HF results to those from 
the harmonic approximation. The deviations are small 
in general and vanish in the strong-dimer limit J' — cx). 
Inside the disordered phase, we observe A'o > 0.7 which 



implies that Uk and Uk < 0.1, see Eqs. (20) and (29 1. The 



fact that -Dk is a small quantity is used in order to derive 
the self-consistent equations (28) and (29), see Appendix 
[C]for details. 

We note that we experienced more difhculties in ob- 
taining self-consistent solutions near the phase boundary 
in the HF approximation as compared to the harmonic 
one, but these problems are cured upon including the 
cubic term, see Sec. |IVB] below. 



B. Hartree-Fock-cubic approximation 

We implement a Hartree-Fock-cubic (HF-cubic) ap- 
proximation by perturbatively adding the effect of 'H3 to 



the mean-field Hamiltonian = Eo + E^" + ^2 + nf^. 



This scheme, which considerably simplifies the treatment 
of the cubic triplet interaction, is similar in spirit to the 
one adopted in Refs. [37] and [39l where an interacting 
Hamiltonian for triplets without the cubic term is treated 
self-consistently first, and then the lowest-order correc- 
tions due to the cubic term are added. 



Using the Bogoliubov transformation ( 14 ) with the 



renormalized coefficients Wk and Uk (27), one can show 



that Ha ^ in terms of the fej^.^ and &kQ operators reads 



X — X — ' 



//V' ^ ^ 



ri(k,p)(6Lpc^l^&k^ + H.c.) 



+ 



k,p 



r2(k,p)(6Lp.^L^-k.+H.c.). (30) 



Here the sum over a,/3,7 has only three components: 
(q;,/3,7) = {x,y,z), {z,x,y), and {y,z,x). The renor- 
malized vertex ri(k, p) is given by 



ri(k,p) 



($k-p - Cp) (Uk-pUpUk + Wk-pWpWk) 

(Ck + Cp) (t^k-pMp^^k + Mk-pt'pMk) 

(Ck-p + ^k) (^^k-pMpUk + Wk-pl^pWk) , (31) 



with being the bare cubic vertex (12), and the ver- 
tex r2(k, p) = — ri(k, p) in addition to the replacements 
Wk ^ Wk. The vertices ri(k, p) and r2(k, p) are illus- 
trated in Fig. |8|i. 

The lowest-order diagrams which contribute to the 
(normal) self-energy I]3(k, a;) are shown in Figs, [sjjb)- 
(e). In each diagram, the solid line corresponds to the 
bare h triplon propagator (we now omit the index a since 
the three triplon branches are degenerate) 



Go-i(k,^) 



ojk + iS. 



(32) 



Note here that no anomalous bare b propagators exist; 
in principle, those are generated in perturbation theory, 
but will be neglected in the following.^" Using standard 
diagrammatic techniques for bosons at zero temperature, 
one shows that only the diagrams (b) and (e) are finite, 



and therefore E3(k, w) = E|''^(k,a;) + T>^'{\i,uj) with 



r?(k,q) 



47V' ^ LJ 
q 



47V 



-E 



ri(k,q) 



t^k-c 



i6 



(33) 



The renormalized triplon energy J7k is given by the 
poles of the full Green's function G(k, w), i.e. 



G-i(k,w) 



Wk - T,3{'k,uj) = 0. 



(34) 



We solve Eq. (34) within the so-called off-shell approx- 



imation, which consists in evaluating the self-energy at 
w — r^k ~ *rk, i.e. 



- iFk " Wk - S3(k, rik - iFk) = 0. 



(35) 



[Recall that in the more common on-shell approximation, 
the self-energy is evaluated at the bare single-particle en- 
ergy, i.e., rik - «fk - <I'k - S3(k,ajk) = 0.] Moreover, 



instead of using Eq. ( 35 ) , causality requires to consider 
f^k-ifk-<;'k-S3(k, rik + ifk) =0. (36) 

The procedure outlined above follows Ref. 35 where 
spin-wave excitations of the isotropic triangular-lattice 
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FIG. 9. (Color online) Triplon dispersion relation along some particular lines of the dimerized Brillouin zone for (a) J' — 3 and 
J" = and (b) J' — 1.5 and J" = 1 at the HF (dashed black line) and HF-cubic (solid red line) approximations. The triplon 
decay rate Fk (dot-dashed blue line) and the bottom of the two-particle continuum (dotted green line) both at the HF-cubic 
approximation are also shown. The A~B line includes the triplon minimum energy of the corresponding approximation. 



AF Heisenberg model (J' — J" = 1) were calculated. 
It is shown by Chernyshev et alW^ that the on-shell 
approximation leads to discontinuities in the spin-wave 
spectrum and concomitant logarithmic singularities in 
the decay rate Fk, and that the ofF-shell approximation 
regularizes such singularities. Finally, the replacement 
f2k — — ^ f^k + in the argument of the self-energy 
guarantees that the quasiparticle pole is in the correct 
half of the complex plane (we refer the reader to the Ap- 
pendix D, Ref. |35t for details). 

In Fig. [9j we compare the renormalized triplon energy 
ilk obtained in the HF-cubic approximation to that from 
HF. Clearly, the inclusion of the cubic term leads to size- 
able changes of the dispersion, which are particularly pro- 
nounced in the incommensurate regime [see Fig.|9]b for 
J' = 1.5 and J" = 1]. For instance, the dispersion along 
r — M is significantly enhanced by the cubic term. Most 
notably, the triplon gap A is renormalized downwards, 
see Fig. [lO] below. In other words, cubic interactions 
tend to destabilize the paramagnetic phase, whereas (re- 
pulsive) quartic interactions have the opposite effect. 

The behavior of the minimum wavevector Qq as a func- 
tion of J" is again qualitatively similar to the harmonic 
one, Fig. [5] Interestingly, the shift of the CIT line due 
to the cubic interaction is opposite (and larger) to that 
induced by the quartic interaction, such that the CIT 
line is now located at some J^it ^.75 which moreover 
depends non-monotonically on J', Fig. [sj 

Fig. [9] also displays the decay rate of triplons caused 
by two-particle decay via "Ha — this is non-zero when- 
ever the single-particle branch is above the two-particle 
threshold. 

We conclude that the cubic triplon term is of vital 
importance for a quantitative description of the triplon 
dynamics in the frustrated coupled-dimer model ([!]), es- 
pecially in the incommensurate regime. This is in con- 



trast to, e.g., the unfrustrated asymmetric bilayer model 
studied in Ref. JiT, where the cubic term leads to minor 
corrections only. A further discussion of physical impli- 
cations of our results is given below. 



V. DISCUSSION 

A. Phase boundary 

For the different levels of approximation, the evolution 
of the energy gap A with the dimerization strength J' is 
shown in Fig. 10 for J" = 0, 0.8, and 1. As expected, 
A increases with J'; the vanishing of the gap defines a 
critical value where the singlet phase becomes unstable 
towards magnetic order. Assuming a continuous QPT, 
we fit the data to 

* ao + aiJ' + a2^ +a3{J'f, (37) 



■J' 

and use the condition A = to estimate the critical cou- 
pling J'^. (Note that the present approximations cannot 
be expected to yield non-trivial critical exponents. Omit- 
ting the 1/ J' fitting term leads to only minor changes of 

J'c- ) 

The resulting phase boundary for < J" < 1.5 is 
shown in Fig. [sj with the HF-cubic result (□) being 
our best approximation. The phase diagram in Fig. [3] 
displays four distinct regions: At large J' we have the 
gapped dimer phase, with spin correlations peaked at 
Qo = (0, 0) for < J" < JciT ^'^'^ incommensurate Qo 
for J" > Jqij-- At the harmonic level J^ix equals 0.75, 
but acquires a J' dependence from anharmonic terms, 
see Figs. [5] and |6] The closing of the triplon gap at 
for a given J" leads to magnetic long-range order 
at wavevector Qo (provided that no other phase inter- 
venes). One concludes that the system displays collinear 




Neel order for J < Jqyt ^'^'^ non-coUinear spiral order for 
J" > JciT- indicated in Fig. [2j an additional param- 
agnetic phase might be realized near the crossing of the 
CIT and the order-disorder transition linespJ with a more 
detailed discussion given in Sec. VB[ 

Let us quantitatively discuss the evolution for two val- 
ues of J" . (i) J" = 0: Here we find a critical coupling 

= 2.65 (2.76) at the harmonic (HF-cubic) level, which 
is quite close to the one determined via QMC simulations 
for the (topologically equivalent) staggered dimerized AF 
Heisenberg model on a square lattice, J'^^^ = 2.5196.'23 
For J' < the syste m ind eed develops collinear Neel 
order with Qq = (0, 0).122IMI (u) j" = where we find 

= 1.09 (1.34) at the harmonic (HF-cubic) level. Again, 
this is a very reasonable result since LRO is certainly 
present at the isotropic point {J' = J" = 1). One differ- 
ence to case (i) is that the ordering wavevector varies with 
J' inside the ordered phase (see appendix [A| for the cor- 
responding classical-limit results); therefore, Qo at the 
phase boundary is distinct from the Goldstone wavevec- 
tors K,K' of the 120° structure, see also Fig. [T|;. 



B. Commensurate— incommensurate transition 



Moreover, in the incommensurate regime near the CIT 
the dispersion along the k2 direction (connecting T and 
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Qo) is anomalously flat. This is illustrated in Fig 
which shows a contour plot of the triplon dispersion fit 
(HF-cubic approximation) for J' = 1.5 and J" = 0.8. 
Qualitatively, the soft dispersion near the CIT will lead 
to an effective dimensional reduction. When the triplon 
gap closes, A = 0, Eq. (38) defines a quantum Lifshitz 
see Sec 



VD below. 



pomt. 

Of course, the physics underlying the CIT is also ob- 
served inside the ordered phase. This has been studied in 
particular in the non-dimerized case, J' = 1. Linear spin- 
wave theory *^ ^ yields a QPT from collinear Neel to spiral 
order at J" = 0.5, where the ordering wavevector changes 
from a commensurate to an incommensurate value. (The 
analysis of the classical ground state for the general case 
J' 1 is given in Appendix|A]) The magnetization curve 
has a minimum at the transition, where the spin-wave ve- 
locity vanishes along one particular direction in k space. 
Coupled-cluster calculation^ find a scenario similar to 
spin- wave theory, but here J" 0.8. Finally, series ex- 
pansion results^ not only indicate a CIT at J" ss 0.75, 
but they also provide evidence for a disordered phase in 
the vicinity of the CIT, i.e, for 0.68 < J" < 0.91. 



As shown above, the minimum wavevector Qo of the 
triplon dispersion in the gapped paramagnetic phase of 
the Heisenberg model ([I]) is locked to (0, 0) for small J", 
while it moves to incommensurate values for larger J" 
(Figs. [5] and [6]), with the boundary being located near 
"^ciT ^0.75. This commensurate- incommensurate tran- 
sition (CIT), driven by increasing magnetic frustration, 
has various consequences. 

Right at the CIT, the quadratic piece of the triplon 
dispersion vanishes in one of the two space directions 
(independent of the level of approximation), i.e., we find 
for the triplon propagator 



G-\k,uj) ^ ~uj^ + + cikl + dik^ 



(38) 



in a small-momentum expansion around the T point, 
where fci 2 are the two components of k perpendicular 
and parallel to the Qo which is taken beyond the CIT. 



C. Two-particle decay 

The cubic triplon term is generically present in 
the mo del (IT] ) and thus enables two-particle decay of 
triplons.'^SEnSuch decay occurs if the one-triplon disper- 
sion moves inside the two-particle continuum. 

As can be seen from Fig. [9] two-particle decay does 
not happen in the commensurate regime of small J", 
where one-particle spectrum is always below the two- 
particle continuum. Note that an exceptional case is 
realized right at the QPT to the magnetically ordered 
state, where the lower bound of the two-particle contin- 
uum coincides with the single-particle branch near the 
ordering wavevector (0,0). This non-trivial coupling has 
been argued to lead to a novel unive rsality class for the 
phase transitioujl^l! see also Sec. VD below. 
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FIG. 11. (Color online) Contour plot of the triplon energy fJk 
(HF-cubic), for J' = 1.5 and J" = 0.8, i.e., near the CIT. 



The situation is different in the incommensurate 
regime, where high-energy triplons can decay into triplon 
pairs, see Fig. [9] for the calculated decay rate. Note that 
such decay never happens near Qo, i.e., cubic interactions 
are not part of the critical theory in the incommensurate 
case. 

As an aside, we note that decay of high-energy modes 
also occurs inside the ordered phase, where spin waves 
can decay into pairs of spin waves provided that the order 
is non-coUinear, for details see Ref. [35l 



D. Quantum criticality 

Let us briefly discuss the quantum phase transitions in 
the phase diagram shown in Fig. [2] 

The transition from the dimer phase to the coUinear 
LRO state with Qo — (0,0) would in principle be ex- 
pected to be in the standard Heisenberg [or 0(3)] uni- 
versality class in Z) = 2 + 1 dimensions; however, the 
cubic triplon term becomes part of the critical theory 
and leads to a new universality class (labelled class B in 
Ref. 



with leading 0(3 ) exp onents and anomalously 

large corrections to scaling P2124| 

The transition from the dimer phase to the non- 
coUinear LRO state does not display this complication 
and therefore is a standard 0(iV) QPT in Z) = 2 -I- 1, 
with = 4 (assuming a single transition to coplanar 
spiral order). 

Finally, we can discuss the multicritical point where 
the collinear LRO, non-coUinear LRO, and dimer phases 
meet. Here, the softness of the triplon dispersion in one 
direction, Eq. ( |38| , implies that the quadratic piece of 
the critical field theory takes the form 



d^xdT^aiS + dl + dl+d^ 



(39) 



where 4>a are the components of the magnetic order pa- 
rameter, and 6 = defines the phase transition point. 
Such (d, m) quantum Lifshitz points (here d = 2, and 
TO = 1 refers to the numb er of "soft" directions) have 
been discussed before,'2S'^but the relevant case of d = 2 



dimensions with undamped order-parameter dynamics 
has not been studied in any detail. A thorough analysis of 
this critical theory is therefore deferred to a future publi- 
cation; here we only make a few qualitative remarks. The 
absence of the quadratic derivative in y direction implies 
that the effective dimensionality is reduced compared to 
a standard 0"* theory, or in other words, the upper crit- 
ical dimension (in the quantum case) is increased from 
d+ = 3tod+ = 3 + to/2 (Ref. While Eq. (IM, sup- 



plemented by a standard quartic term (note that a cubic 
ternPll rnight be important as well), describes a continu- 
ous multicritical point, it is conceivable that this is pre- 
empted by a fluctuation-induced transition into a novel 
phase. We speculate that this could be a non-trivial para- 
magnetic phase (e.g. with a symmetry-breaking dimer- 
ization), as deduced for J' = 1 from series-expansion 
studies.!^ 



E. Application to EtMe3P[Pd(dmit)2]2 

As mentioned in the Introduction, the experimental 
finding^^mil] indicate that the organic Mott insulator 
EtMe3P[Pd(dmit)2]2 realizes a columnar VBS phase at 
low temperature and pressure. According to Ref. [22l 
this system can be described by the Heisenberg model 
^ with J' = 1 and J" = 1.05. Although this con- 
figuration is outside of the VBS region predicted by our 
bond-operator analysis, it is quite close to the VBS phase 
boundary (see the orange diamond. Fig. [3]). It is conceiv- 
able that a combination of longer-range or ring-exchange 
interactions, which arise in the proximity to the Mott 
transition,!^ increase the level of magnetic frustration, 
which is then released by a lattice dimerization through 
magnetoelastic couplings (note that organic compounds 
of the X[Pd(dmit)2]2 display a rather soft lattice). As 
a result, a VBS phase can emerge, whose magnetic cou- 
plings are explicitly dimerized as in our Hamiltonian ([T]) . 

Therefore we believe that the excitation spectrum of 
the organic compound may display features qualitatively 
similar to the ones shown, e.g., in Fig. |9]b, which corre- 
sponds to a configuration (J' = 1.5, J" = 1) close to the 
critical line. Hence we predict the spin correlations in the 
VBS phase of EtMe3P[Pd(dmit)2]2 to be incommensu- 
rate, which may be checked in future neutron scattering 
or NMR experiments. Moreover it would be interesting 
to see whether one could drive the material into a state 
with magnetic LRO, by moving farther away from the 
Mott transition (and thus reducing the influence of mag- 
netic couplings beyond nearest-neighbor exchange). 



VI. SUMMARY 

We have studied a dimerized Heisenberg AF on a spa- 
tially anisotropic triangular lattice, with focus on its 
quantum paramagnetic (i.e. dimer) phase. Starting from 
bond-operator mean-field theory, we have included inter- 
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action corrections to the harmonic approximation at the 
Hartree-Fock level for the quartic term and in second- 
order perturbation theory for the cubic term. We have 
shown that this Hartree-Fock-cubic approximation gives 
sensible results for the triplon dispersion for the investi- 
gated part of the parameter space (away from the one- 
dimensional limit). The resulting boundary of the dimer 
phase where the triplon gap closes, Fig. [3j is one of our 
main results; its location in the unfrustrated limit J" = 
is in good quantitative agreement with QMC data. 

The minimum wavevector Qg of the triplon spectrum 
displays a commensurate-incommensurate transition as 
frustration is increased. At this CIT the quadratic piece 
of the triplon dispersion vanishes in one momentum direc- 
tion. The closing of the triplon gap at the CIT leads to a 
quantum Lifshitz point, the detailed study of which is left 
for future work: it either leads to distinct power-law crit- 
ical behavior or it may even be inherently unstable, lead- 
ing to a novel intermediate phase, Fig. [2j Remarkably, 
we find the critical J' to be only slightly above unity, sug- 
gesting that the origin of the non-trivial VBS phase seen 
in series-expansion studies^ for the non-dimerized model 
is indeed the enhanced fluctuations near the quantum 
Lifshitz point. Finally, in the incommensurate regime 
of the dimer phase the cubic triplon interaction leads to 
two-particle decay of triplons at high energies. 

Our study paves the way to further investigations of 
frustrated dimerized magnets. As those are difficult to 
access using QMC simulations, due to the inherent minus 
sign problem, bond-operator as well as series-expansion 
techniques are often the methods of choice. Clearly, the 
inclusion of longer-range and cyclic exchange interactions 
on triangular and Kagome lattice geometries would be 
most interesting. On the methodological side, a self- 
consistent treatment of the cubic term might further 
improve the numerical accuracy. In addition, the for- 
mation of bound states of triplons - which would sig- 
nify the instability towards phases other than those with 
magnetic LRO - should be studied; we leave this for fu- 
ture work. Finally, the interplay of magnetoelastic cou- 
plings and longer-range exchange interactions should be 
studied near the isotropic point, with an eye towards 
EtMe3P[Pd(dmit)2]2- 
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Appendix A: Classical phase diagram 

To determine the classical phases of the model ([!]), we 
parametrize the spins S,^ and (Fig. [l^) according to 

S,J = ei cos(Q • Ri) + 62 sin(Q • R^), 

(Al) 

— 63 cos(Q • Ri) -f 64 sin(Q • R^), 



which assumes coplanar order. Here e, 
vectors which obey: ei • 62 = 63 • 64 = 
cos9, 62 • 63 = — ei • 64 = sinfl. Substituting Eq. (Al) 



are a set of unit 
ei -63 = 62 -64 = 



into the Hamiltonian Q, it is easy to see that the total 
energy is E = {N/2)S^q{9) with 

Jq (9) = J' cos 9 + 2 J" cos{Q^) + cos(Q^ + 9) 

+ cosiVSQy + 9)+ cos(g^ + VSQy + 9). (A2) 

For fixed J' and J", the ground state energy is deter- 
mined by minimizing Eq. ( A2 ) with respect to the com- 



ponents of the vector Q and the angle 9. 

a) collinear order, with 



We find two phases. Fig. 12 



a commensurate ordering wavevector Q = (0, 0) and 
9 = TT, realized for J" < J'cn^iJ'), and b) non-collinear 
order, with incommensurate Q — {Qxm,Qym) and 9 = 
71" ~ Qxm/2 — \/3Qym- The components of Q solve 

^ sin{Qx/2) cos{ViQy) - sin{Qx/2 - ViQy) = 0, 



4 J" cos{Q^/2) ~ ncos(V3Qy) = 1, 



(A3) 



with /i = 2J'/(1 + J'). Numerical solutions of (A3) can 



be easily obtained, which yield the phase boundary in 
Fig.[l2| 

Interestingly, it is possible to solve Eqs. (A3) analyt- 
ically for the particular cases J' ==0, 1, and od which, 
respectively, correspond to /i = 0, 1, and 2. The solutions 
can be written in the following way: 

, ,-,.2 

2 arccos 



l-fl'^+2^l 
4J" - m(m - 1) 



(A4) 



/ l-/i2 + 2^ 
— ^ arccos r 

V3 



Here, Aq — — A2 — 1 and Ai = 0. Since the argument of 
arccos(x) is within the range — 1 < a; < 1, we conclude 
that the non-collinear phase is stable for J" > Jqyt where 



T" 



J' 

2(1 + J')' 



(A5) 



valid for J' = 0, 1, and 00. (For J' = 00, we obtain 



T" 



CIT 



0.75 and 



i.e., the classical re- 



sult matches the corresponding features of bond-operator 
result for S = 1/2.) 

Quantum corrections to the classical phase boundary 
can be obtained using non-linear spin-wave theory. For 
J' — 1 this has recently been demonstrated to move the 
CIT boundary to J" = 0.77 at order l/S (Ref.dli)- 
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FIG. 12. (Color online) Classical phase diagram for the 
Heisenberg model ([T]), obtained from numerically solving 
Eq. (TASI). 



Appendix B: Effective triplet Hamiltonian in real 
space 

Here we quote the full expressions for the real-space 
Hamiltonian ([6|, which can be obtained using straight- 
forward algebra via Eq. ([3|. 



7^0 



3J' 



E 



(Bl) 



J' 



('^j'^j+n^i a^'J-l-n a + H.C.) 



(B2) 



-HH.c.]- (joj+n)}, (B3) 

where the runs over the sites of the dimerized lattice, 
and the gmin) coefhcients read 

54(«) = (2J" l)5„,i + (5„,2 + (Jn.a- 
Appendix C: Self-consistent equations 



Here we derive the self-consistent Eqs. (19)-(20l in- 
volved in the harmonic approximation, and Eqs. (28|- 
(291 related to the HF one. As mentioned in Sec. HI 



the starting point is the saddle point conditions: (i) 
dE^/dNo = and (ii) ^Ea/^^JL = 0. 



Let us consider the harmonic case. From the expres- 

;y, 
9w„ 



sion ( 16 ) for the ground state energy, we have 
5^0 



q 



(CI) 



Since dA^/dNo = dBjdNo = BJNq, see Eqs. ^ 
(11), condition (i) yields Eq. (19). Moreover, 



\n{n,-i) + Iy. 



(C2) 



One can easily see that dAq/dfi = — 1 and dcuq/d I-l = 
— Aq/wq, and therefore condition (ii) leads to Eq. (20). 



Turning to the Hartree-Fock approximation, the equiv- 
alent of Eq. ( CI ) is here given by 



q 



dNo 







dAq 



(C3) 



The contribution due to E^^ , Eq. (l24 



since it is ©(Wq) and therefore smal 



other terms (see discussion in Sec. IV A). The last term 
of the above equation can be written as 



, can be neglected 
compared to the 



/ = 



2 ^ 
q 



-1 



A„\ dA, 



Sq dBq 



Wq J dNo Wq dNo 



with 



dAq 

dNo 

dB„ 



Sq 

No 
B„ 



7k-q 



N' ^ 

k 

]^ E "^k-q 



d 

d 
dN~o 



Bk 



dNo No 
Assuming that 

dAJdNo - dBjdNo « BjNo 
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and using the expressions (27), we show that / ~ 
(3/A^o)Eq^q^q(^^q-'«q). and therefore Eq. ([28| fol- 
lows from condition (i). Some words about the assump- 
tion (C4) are here in order: we realize that it is a good 



approximation as long as Vq is small. Indeed, using the 
relations (27), it is possible to show that 
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Substituting ( C4 ) in the above expression, one can easily 
see that dA^/dNo ~ Bq/No + 0{v\s). Similar expression 
holds for dBq/dNo. 



The analog of Eq. ( C2 ) reads 



(C5) 
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Here 
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In order to be consistent with the above approximation, 
we assume that dAq/dfi w — 1 and dBq/dfi ~ which 
imphes duj^/d^ « — Aq/wq. Notice that this is indeed 
a reasonable approximation because condition (ii) yields 
to nothing else but the conservation (on average) of the 
total number of bosons per site of the dimerized lattice, 
i.e., Eq. 
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